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Spin-Orbit and Tensor Forces in Heavy-quark Light-quark Mesons: Implications of 

the New State at 2.32 GeV 
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We consider the spectroscopy of heavy-quark light-quark mesons with a simple model based on the 
non-relativistic reduction of vector and scalar exchange between fermions. Four forces are induced: 
the spin-orbit forces on the light and heavy quark spins, the tensor force, and a spin-spin force. If 
the vector force is Coulombic, the spin-spin force is a contact interaction, and the tensor force and 
spin-orbit force on the heavy quark to order l/m 1 m2 are directly proportional. As a result, just 
two independent parameters characterize these perturbations. The measurement of the masses of 
three p-wave states suffices to predict the mass of the fourth. This technique is applied to the D s 
system, where the newly discovered state at 2.32 GeV provides the third measured level, and to the 
D system. The mixing of the two J p = 1 + p-wave states is reflected in their widths and provides 
additional constraints. The resulting picture is at odds with previous expectations and raises new 
puzzles. 

PACS numbers: 12.39, 13.20Fc 



I. INTRODUCTION 

Mesons composed of one light quark and one heavy 
quark are quite analogous to the hydrogen atom [1] and 
can be analyzed using the traditional methods [2]. The 
fine and hyperfine structures have direct analogues in 
meson spectroscopy, but the confinement of quarks re- 
quires that potential cannot be purely Coulombic. A con- 
venient phenomenological approach is to postulate that 
there are two separate static potentials, one arising as the 
zcroth component of a vector potential and the other as 
a Lorentz scalar. Asymptotic freedom suggests that the 
vector potential be Coulombic, while confinement sug- 
gests that the scalar be a linear potential. Such mod- 
els have had reasonable success in describing the spec- 
troscopy of the D, D s , B, and B s systems [3-7]. 

The recently discovered state with a mass of 2.32 GeV 
decaying to D s ir a [9] appears to be a cs p-wave meson 
with J p = + . Typical predictions for its mass were 
near 2.5 GeV, above the threshold for the strong decay to 
DK. Here we analyze the p-wave cs and cu/d mesons in a 
simple model that is more general, though less predictive, 
than the models described above. It is more general in 
that we do not specialize to a scalar potential that is 
linear. 

Our concern here is restricted to the p-wave states. 
While our primary interest is in the D s system, we con- 
sider as well the analogous non-strange D mesons. It is 
not only the mass spectrum that needs to be addressed, 
but also the pattern of decay widths. The early studies, 
which predicted a much higher mass for the J p = + cs, 
were successful in explaining the narrow width of the ob- 
served J p = 1 + states in both the D and D s systems. We 



find that the combined constraints of the mass spectrum 
as now known (though lacunae remain) and the decay 
patterns make it difficult to find a consistent picture of 
the p-wave charmed strange and non-strange mesons. 

By analogy with the hydrogen atom, a convenient clas- 
sification of states is given in terms of j = £ + s±, where 
£ is the orbital angular momentum and Si is the spin of 
the light quark. The light-quark angular momentum j 
is conserved in the limit in which the heavy-quark mass, 
J7i2, goes to infinity. For p-wave states, the values of j 
are 3/2 and 1/2. These levels are split by the ordinary 
spin-orbit force. The hyperfine structure is proportional 
to \jm,\m2 and includes the spin-orbit coupling of the 
heavy quark, a spin-spin interaction, and the tensor force. 
All these contribute to further reducing the degeneracy 
so that J = j + s 2 is conserved, but not j alone. The 
result is four distinct states with J p = + , 1 + , 1 + , and 
2 + . The two J p = 1+ states are mixtures of j = 3/2 and 
j = 1/2. The J = 2 and J = states are pure S = 1 
states, where S = si + S2 is the total spin. The J = 1 
eigenstates are mixtures of spin-triplet and spin-singlet. 

The p-wave D mesons decay by pion emission. The 
2 + state decays through d-wave emission to the ground 
state _ , D or to the 1 _ , D* . The 1 + states can decay, in 
principle, either by d-wave or s-wave pion emission. The 
limited phases space favors the s-wave decay. However, 
the decay of j = 3/2 state to the D* , which has j = 1/2, 
is forbidden to the extent that j is conserved. Indeed, 
the J p = 1+ state at 2.422 GeV is narrow. 

For p-wave D s mesons, decay by pion emission is for- 
bidden by isospin invariance. The D s states at 2.572 GeV 
and 2.536 GeV decay to DK and D*K. In a fashion anal- 
ogous to the pattern in the D system, the J p = 1 + state 
at 2.536 GeV is also narrow. 
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II. SPECTRUM FOR P STATES 

We base our analysis on the quasi-static potential, in- 
cluding the spin-dependent forces, which computed from 
the Lorentz-invariant fermion-antifcrmion scattering am- 



plitude using Feynman diagrams, replacing the vector or 
scalar propagators by the Fourier transforms of the pos- 
tulated potentials, S(r) and V(r). In this way, neglecting 
velocity-dependent, but spin-independent terms, we find 



Am\ Am\) \ r J \ 2m\rn 2 
+ (— - V") S 12 + t-^—^V tn-v* (1) 



where the tensor force operator is 

S\2 = 3<ri-r cr 2 r — a\-cr 2 . 



(2) 



We imagine solving the corresponding differential 
equation with the potential V + S, then computing 
the fine structure and hyperfine structure perturbativcly. 
We work consistently only to order \jm 2 . In general, 
there are four independent matrix elements to consider, 
((V - S')/r), (V'/r), (-V" + V'/r), and (X7 2 V). How- 
ever, if V is Coulombic, -V" + V'/r = W'/r and (VV) 
vanishes except for s waves, where it gives a contact in- 
teraction. It follows that for Coulombic V and for I > 0, 
there are only two independent matrix elements. Three 
measured p-wave masses will give two splittings, which 
will determine the matrix elements and allow us to pre- 
dict the mass of the fourth p-wave state. 

The mass operator for the splittings of the four states 
of any multiplet of orbital angular momentum different 
from zero can be written (to order l/m 2 ) [8] 



M = \£ ■ si + 4r£ • s 2 + tS 



with the notation 
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In practice, we shall apply this operator to p-wave states 
for the D or D s system. Henceforth, we will use A and 
r to indicate the expectation values of the quantities in 
Eq. (4). The values of A and r will be different for the 
D and D s systems. For the assumed attractive, Coulom- 
bic V, the tensor-force energy r is manifestly positive. 
However, the spin-orbit energy A can be either positive 
or negative, depending on the relation between the po- 
tentials V and S. 

Because there is no single basis that diagonalizes all 
the interactions, we need to fix one basis, then calculate 
the mixing between the two J = 1 states. We choose the 



basis in which j 2 is diagonal. The eigenstates | J, j, m) 
of J 2 , j 2 and J z can be written in terms of eigenstates of 
J 2 , S 2 , and J z , or in terms of eigenstates of J 2 , (j 1 ) 2 = 
(£ + s 2 ) 2 and J z . The details are given in the Appendix. 

Up to an additive constant common to the four p-wave 
states, the masses of the J = 2 and J = states are 



A 8 
M 2 = - + -t 



M = —A — 8r 



(5) 
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while the masses of the two J = 1 state are obtained by 
diagonalizing the matrix in the | J, j, m) basis 



/ a _ 8 

'2 3 ' 



-A+|r 



The two eigenmasses for J = 1 arc then 



Mi± 



A 



± 



(6) 



(7) 



The eigenmasses are shown as functions of A/t in Fig. 1. 
Also shown is the j = 1/2 fraction of the higher mass 
J p = 1+ state. The vertical bands correspond to the- 
oretical models and to data, as explained in the next 
section. 

If we define the mass splittings among three of the four 
states, 



D 2 = M 2 - M 
Di = Mi - M 



(8) 



we find 



If) 2 

— D 2 D 1 

87 2 29 1 



± V(^-| D ') +55<"I-W»> 

»-|ft-fr (9) 

A very strong scalar potential S leads to "inversion," 
namely the j = 1/2 states lying above the j = 3/2 states. 
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III. APPLICATION TO THE D AND D 3 

SYSTEMS 

Using the formulae of the preceding Section, we can 
use the measured masses of three p-wave states to predict 
the fourth. Alternatively, we can take three masses that 
are predicted theoretically and confirm that we find the 
fourth predicted mass. Applying this to the D and D s 
systems as given in Ref. [7] we indeed find congruence 
and determine the values of A and r shown in Table I. 
The leftmost vertical bar in Fig. 1 indicates the range 
found in Ref. [7] . The negative value of A is indicative of 
inversion. 

The result of inversion is that the higher-mass J p = 1 + 
state actually lies above the 2 + state. This prediction is 
contradicted by the reports from Belle [11] of J p = 1+ 
state at 2.400 GeV. Using the well established D states at 
2.459 GeV and 2.422 GeV, and the reported J p = 0+ at 
2.290 GeV in our Eq. (3) leads to two solutions, labeled 
A and B in Table I. Both give masses near the observed 
state at 2.400 GeV. From the lower graph in Fig. 1 wc 
see that the solution with the higher value of A/r results 
in the higher-mass J p = 1 + state (2.422 GeV) being 
nearly entirely j — 3/2, while the lower value (Solution 
A), would give it nearly equal contributions from j = 1/2 
and j — 3/2. 

We can differentiate between the two solutions, both 
for D and D s by considering the widths. The exper- 
imental widths and theoretical estimates are shown in 
Table II. The theoretical estimates are obtained from 
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FIG. 1: Upper: the energy-levels for the four p-wave states 
as a function of the ratio of the spin-orbit to tensor energies, 
in units of the tensor energy. Lower: the percentage of the 
more massive J F = 1 + state that comes from the j = 1/2 
state. The D masses are from [10, 11] and the D s masses 
from [9, 101. 



TABLE I: Masses of the various p-wave states in the D and 
D a systems and the spin-orbit and tensor energies. The ex- 
perimental masses for the D states at 2.400 GeV and 2.290 
GeV are from Ref. [11]. The mass of the D s state at 2.317 
GeV is from Ref. [9]. The values of A and r for Ref. [7] were 
obtained by fitting their mass spectrum with the ansatz of 
Eq. 3. The square brackets indicate values that were used as 
inputs to the fits that determined the remaining mass and the 
values of A and r. At the time of Ref. [7] only the masses of 
the 2 + states and the Di(2420) and D s i(2535) were known. 





Exp. 


Theory 




Ref. 


[9-11] Sol. A Sol. B Ref. [7] 


D mesons 


M(2+)(GeV) 


2.459 [2.459] 


[2.459] 


2.460 


M(l+)(GeV) 


2.400 2.400 


2.385 


2.490 


M(l+)(GeV) 


2.422 [2.422] 


[2.422] 


2.417 


M(0+)(GeV) 


2.290 [2.290] 


[2.290] 


2.377 


A (MeV) 


39 


54 


-11 


t (MeV) 


11 


9 


11 


D s mesons 


M(2+)(GeV) 


2.572 [2.572] [2.572] 


2.581 


M(l+)(GeV) 


2.480 


2.408 


2.605 


M(l+)(GeV) 


2.536 [2.536] 


[2.536] 


2.535 


M(0+)(GeV) 


2.317 [2.317] 


[2.317] 


2.487 


A (MeV) 


43 


115 


-7 


r (MeV) 


20 


9 


11 



TABLE II: Decay widths of p-wave D and D a states in 
MeV. The theoretical values are derived from [7] using phase- 
space corrections to adjust for the masses known now. The 
widths shown for the 1)2(2460) states are obtained from 
the total width of 23 ± 5 0~ and the measured [10] ratio 
r(L> + 7r~)/r(D*+7r~) = 2.3±0.6. For the D sJ (2573) we have 
assigned the entire width to the decay to D(1865)A' since the 
decay to D* (2007)K has not been seen. 







Exp. 


Theory: 






[10, 11] 


s-wave d-wave 


D mesons 


L>a(2460) -> 


L>(1865)tt 


16 ±4 


16 


L^ (2460) -> 


D*(2007)tt 


7±3 


9 


Di(2422) -> 


D*(2007)7r 


18.9+ts 


90 10 


Di(2400) -> 


D*(2007)7r 


380 ± 100 ± 100 


100 


L»S(2290) 




305 ± 30 ± 25 


100 


D a mesons 


L>a(2573) -> 


D(1865)_K: 


15 1$ 


9 


L»a(2573) -> 


D*(2007)K 




1.4 


Di(2535) -> 


D* (2007) A 


< 2.3 


100 0.3 



[7] after making phase-space corrections. The widths for 
s-wave and d-wave decay are shown separately for the 
J p = 1+ states. The proper combination depends on 
the mixing of the j = 3/2 and j = 1/2 states. 

Referring to Fig. 1 and Table I, we see that the model 
of [7], which has A < 0, has the J = 1 state that is nearly 
pure j = 1/2 above the J — 2 state. The lower J = 1 
state is identified with the £>i(2420). Because it is nearly 
entirely j — 3/2 its decay by pion emission must be d- 
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wave. This causes a suppression that conforms to the 
measured width. However, the mass spectrum predicted 
by Ref. [7] is not in good agreement with the data. Of 
our two solutions for the D data, only one has a large 
value of A/r as needed to make the lower mass J = 1 
state broad, as required by the data. 

Turning to the D s system, we find a similar situation. 
While Ref. [7] would have A < 0, the known masses, 
which now include the + state at 2.32 GeV, require 
A > 0. Indeed, to suppress the width of the D s (2535), a 
large A is needed. This is the case for Solution A, where 
the mass predicted for the fourth p-wave state is 2408 
MeV. One of us ( JD J) notes that in Ref. [9] , there is an 
apparent signal in the invariant mass of D s tt q ^ 2.46 GeV. 
This would be consistent with a J p = 1+ state decaying 
through L>* 7 (2317)7 or through L>*(2112)tt . The mass 
fits better with Solution A of Table I, but the narrow 
width of the £> s i(2535) favors Solution B. 

Both the D and D s systems show deviations from the 
pattern anticipated by potential models. While "inver- 
sion", i.e. A < 0, has been a favored prediction, it is not 



in agreement with the data. This suggests that the ansatz 
taken for the potentials V and S may not be as simple as 
assumed. We have allowed S to be any potential, while 
requiring that V be Coulombic; Ref. [7] required that 
S be linear. Relaxing that restriction might lead to im- 
proved agreement with the data for this highly predictive 
approach, which actually uses the solutions to the Dirac 
equation. On the other hand, even with the flexibility 
allowed in our approach, the resulting picture is not en- 
tirely attractive. We expect the confining potential to be 
quite important since the p-wave states are not concen- 
trated at the origin. This would lead to a suppression of 
A through the contribution of the —S'/r term. This is 
not borne out by the values of A we deduce. 

The discovery of the D*j(2317) has provided an im- 
portant clue to heavy-quark light-quark spectroscopy by 
nailing down a p-wave state with j = 1/2. Puzzles re- 
main. The anticipated discovery of the accompanying 
j = 1/2 state with J = 1 should add important new in- 
formation, but it is not likely to resolve all the questions 
we have described. 
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IV. APPENDIX 

Ions betweeen the three bases that diagonalize j 2 , j' 2 , and S 2 are 

= \J = £+1, S = 1, to) 
J = £, j = £-1/2, to) = yfg±\ J = £,S = l,m)- yjrf^l J = £,S = 0,m) 

= 1, m) + J = £,S = Q,m) 



J = l,j = £+1/2, to) =^Z\J = £,S=l, 
J = £-l, j = £-l/2, to) = \J = £- 1, S = 1, to) 



= £+1, j' = £+1/2, to) = \ J = £+ 1, S= 1, to) 

J = £, f = £ - 1/2, to) = | J = £,S = l,m) + yfJ^\J = £,S = 0,m) 

J = £,f = £+ 1/2, to) = y /^\J = £ t S=l,m)- y fg±\J = e,S = 0,m) 
J = £-1,3' = £-1/2, to) = \ J = £- 1, S= 1, to) 

J = £+l,j = £+1/2, to) = \ J = £+ 1, j' = £+1/2, to) 

J = 1,3=1- 1/2, m) =^ ri \J = £,j' = £- 1/2, to) + V 2 j ( ^ +1) |J = £, f = £ + 1/2, to) 

J = £,3=£+ 1/2, to) = ^gjP | J = £, 3' = t - 1/2, m)- J j ri \J = £, 3' = £ + 1/2, to) 
J = £-l,j = £-1/2, to) = \ J = £- 1, j' = £-1/2, to) 

In these three bases it is easy to evaluate the matrix elements of t-s\, £-S2 

(Jjm \2£- Sl I Jjto) = 3(3 + 1) - £(£ + 1) - 3/4 
(Jj'm\ 2l-s 2 \Jj'm) = j'(j> + 1) - £(£ + 1) - 3/4 



and of S12 (all of whose matrix elements vanish if the initial or final state has S = 0) 
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(J = £ - 1, 5 = 1, m|5i 2 | J = ^ - 1, 5 = 1, m) = -2^r 
(J = £,S = l,m \S 12 \J = £,S =l,m ) = +2 



(J = ^ + 1,5= l,m|5i 2 |J = £ + 1,5 = l,m) 



J 2,7+1 



from which wc find 



and 



(J = 


e 


+ l,j=£ + l/2,m 


2i 


S2 


J = £ + 


i,i 


= £+l/2,m) 


(J= 




j =£+l/2,m 


2£ 


S2 


J= 


1,3 


= £+l/2,m) 


(J= 


e, 


j = £-1/2,171 


2£ 


S2 


J = 


t,3 


= £-l/2,m) 


(J= 


e, 


j=£-l/2,m 


2£ 


S2 


J = 


i,3 


= £+l/2,m) 


(J = 


£ 


-l,j=£-l/2,m 


2£ 


S2 


J = £- 


1,3 


= £-l/2,m) 



(J = 


£ + 1, 


3 = 


£+1/2, m 


5l2 


J = £ + 


1,3 


= £+l/2,m) 


(J = 


t,3 = 


£ + 


1/2, m 


5l2 


J = 




= £+l/2,m) 


(J = 


e,j = 


£- 


1/2, m 


5l2 


J = 


e,j 


= £-l/2,m) 


(J = 


1,3 = 


£- 


1/2, m 


5l2 


J = 


*,3 


= £+l/2,m) 


(J = 


£-1, 


i = 


£-l/2,m 


5l2 


J = £- 


1,3 


= £-l/2,m) 



J -I 

J(2J+3) 
2,7+1 

(2J-1)(J+1) 
2,7+1 



2^/.7(./+l) 
2J+1 

-J-2 



J 2,7+1 



— 2,7 



2J+1 

2(.7+l) 
2J+1 

2^/./(J+l) 
2J+1 

O J+2 
Z 2,7+1 
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